function [  ] = splitStep(  )
   close all
    
        
    % initialize simulation parameters
    dimR = 200;
    dimZ = 2001;
    z= 0.5;
    h= 1e-9;
    dt = 1e-15;
    
    % set pulse
    pulseDuration = 50e-15;
    a0 = 1e20;
    dr = 0.5e-6;
    w0 = 20e-6;
    k = 1e26;
    u = createPulse (dimR,dimZ, a0, dr, w0,dt,pulseDuration);
    
    % set material coefficents
    alpha =  0; 1e7;
    b2 =   1e-19; %2nd order disp. (s2/m)
    gamma= 1e-31;% non linearity in /W/m
    
    
    % initialize spectral steps
    dw=1/dimZ/dt*2*pi;
    w=(-1*dimZ/2:1:dimZ/2-1)*dw;
    w=fftshift(w);
    
   % main loop  
    for jj=h:h:z
        jj
        
        
        u = splitStep1h( u , dimR, h ,w, alpha, b2 , gamma,k,dr);
        
        %figure (1);
        %wl = [1:1:2001]-1000;
        %plot (wl, fftshift (abs(fft (u(50,:)))));
       
        % plot spati-temporal domain
        figure(1);
        colormap hot;
        imagesc (abs(u));
        colorbar('SouthOutside');
        drawnow
        
    end
    
 
   
end
    





